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Abstract 

Diffusive transport is a universal phenomenon, throughout both biological and physical sciences, and 
models of diffusion are routinely used to interrogate diffusion-driven processes. However, most models 
neglect to take into account the role of volume exclusion, which can significantly alter diffusive transport, 
particularly within biological systems where the diffusing particles might occupy a significant fraction of the 
available space. In this work we use a random walk approach to provide a means to reconcile models that 
incorporate crowding effects on different spatial scales. Our work demonstrates that coarse-grained models 
incorporating simplified descriptions of excluded volume can be used in many circumstances, but that care 
must be taken in pushing the coarse-graining process too far. 
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I. INTRODUCTION 


Throughout the physical and biological sciences, diffusive transport is ubiquitous, and it takes 
places on a wide range of spatial and temporal scales. For example, in biology diffusion is a key 
transport process that regulates events and interactions on levels ranging from those describing the 
behaviours of ions and subcellular macromolecules, to those of cells, tissues and organisms [T]. Less 
well understood, however, is how volume-exclusion-driven crowding impacts upon these diffusive 
transport processes, despite the inherent fact that all diffusing particles exclude other particles 
from occupying the same region in space mi- 

The majority of models of diffusive processes neglect to take into account excluded volume ef¬ 
fects. The predominant models describing diffusive transport over a range of spatial scales, and with 
varying excluded volume fractions, are ‘diffusion’ partial differential equation (PDF) models with 
a constant diffusion coefficient [IHS]) and random-walk-based models of point particles [71[H], both 
of which entirely neglect the impact of volume exclusion. Other models include phenomenological 
descriptions of volume exclusion effects by imposing that, for example, the diffusion coefficients 
of PDF models or the ‘jump rates’ associated with random walk models depend locally on the 
particle density [MS]. However, these phenomenological descriptions are usually chosen on an ad 
hoc basis, and the ramifications of choosing one description over another rarely explored in detail. 

In this work we employ the framework of a lattice-based random walk to explore how volume 
exclusion may be taken into account at different spatial scales, and describe how to define the 
jump rates of the random walkers so as to provide a consistent description of the effects of volume 
exclusion across spatial scales. An additional advantage of our approach is that it provides for 
significant time savings in the computational simulation of volume excluding individual-particle- 
level models of diffusive transport. 


II. MODELLING DIFFUSIVE TRANSPORT ON DIFFERENT SCALES 

For simplicity, we consider a one-dimensional, lattice-based random walk model of diffusive 
transport in which the motile particles have length h. We work with the domain x G [0, L], where 
L = Nh for some V G N, so that the domain can hold at most N particles. We impose a uniform 
lattice consisting of K compartments on the domain, so that the length of each compartment is 
L/K, and we work only with choices of K for which m = Lj{Kh) is a positive integer. This means 
that N = Km and at most m particles fit into each of the K compartments UMM- We model 
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diffusion as a series of jumps between compartments, and impose zero-flux boundary conditions on 
the domain (so that particle number is conserved). 
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FIG. 1. Representation of particle positions within compartments of different capacity. Shaded cells rep¬ 
resent particles, and white cells represent unused capacity. As m increases, the spatial resolution coarsens, 
and this work reconciles the descriptions at these different scales. 


We move between different levels of spatial resolution by varying the compartment capacity, m: 
smaller values of m resolve changes in particle density on a finer spatial scale than larger values. 
The two limiting cases are: full exclusion, m = 1, so that compartments contain at most one 
particle (and the position of the particle is fixed); and no exclusion, m —>■ oo, so there is no limit on 
the number of particles per compartment. In this work, we shall term the m = 1 case ‘accurate’, 
in the sense that no assumptions are made on the positions of particles with each compartment. 

The common, phenomenological approach taken in the literature is to define the jump rates 
between compartments as 
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where is the number of particles in compartment j when each compartment has capacity 
m m- The function /I™) describes the effects of volume exclusion, effectively specifying the 
proportion of jumps that ‘fail’ due to crowding [18]. The scaling of the jump rate with the square 
of the compartment size {mh) can be justified from mean first passage time approaches |19] . 
Effects such as adhesion are often included by assuming the Tf to also be a function of and 
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Typically, these models on different scales are interrogated either by: (i) using repeated simula¬ 
tion of the random walk models to estimate summary statistics of interest; (ii) deriving and solving 
ordinary differential equations (ODEs) for the expected particle number per compartment; or (iii) 
deriving PDE models in the limit /i —)• 0 and using standard analytical and numerical techniques 
for PDEs to explore model behaviours. However, to the best of our knowledge, there has been 
little exploration of the effects of choosing different functional forms for on various summary 
statistics of the random walk models as the compartment capacity, m, varies. As such, one of the 
aims of this work is to understand how the mean and variance of particle numbers changes as we 
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move across spatial scales (by varying m) and to provide a systematic derivation of coarse-grained 
{m > 1) models from the accurate (m = 1) model. 

Sensible choices of require: 


1. the volume exclusion function to be zero when the compartment is empty, = 0; 

2. the volume exclusion function to be unity when the compartment is at capacity, (m) = 1. 


One of the goals of this work is to elucidate functional forms for /(’"I that give rise to behaviours 
that are conserved across spatial scales. We do this by considering equations for mean and variance 
of particle numbers. 

Mean particle numbers. The evolution of mean particle number in the compartment 
when the transition probabilities are as in Eq. 0 is given by the ordinary differential equation 
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for 2 < j < K — 1, where (•) denotes the expectation and Similar expressions 

apply for the boundary compartments, j = \,K. 

The case m = 1. Here at most one particle can occupy each compartment, the conditions 
stated above are enough to fully define and we have m 
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for 2 < j < K — 1. Similar expressions apply for the boundary compartments, j = l,K. Eq. 0 
is a semi-discrete diffusion equation and, in the limit /i —)• 0, it gives rise to the diffusion equation 
with constant diffusion coefficient, D. 

The case m = oo. Letting m —)• oo entails the limit of zero volume exclusion. To analyse the 
evolution of mean particle number we represent the compartment size as h = L/K so that, similar 
to the m = 1, case we have 
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for 2 < j < K — 1. Similar expressions apply for the boundary compartments, j = 1,K. As for 
the m = 1 case, Eq. 0 is a semi-discrete diffusion equation and, in the limit h —0, it gives rise 
to the diffusion equation with constant diffusion coefficient, D. 
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The case 1 < m < oo. To ensure the model is consistent across spatial scales, it is appropriate 
to confine choices for for 1 < m < oo to those that also give rise to a semi-discrete diffusion 
equation with constant diffusion coefficient for mean particle numbers. The only choice for /I™) is 
then 


which gives, as anticipated. 
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dt rn^K^ 

for 2 < j < K — 1. Similar expressions apply for the boundary compartments, j = \,K. 

Variance of particle numbers. For the choice of volume exclusion function given in Eq. © 
we can also obtain equations for the evolution of the variance of particle numbers: 
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for 2 < j < K — 1, where is the variance of particle numbers in compartment j, and is 
the covariance of particle numbers in compartments j and k: 
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for l<j<k — 1<K, l<k + l<j<K. 
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(9) 


Similar expressions can be found for the boundary compartments, j = l,K. 


III. CONSISTENCY OF THE CHOICE OF VOLUME EXCLUSION FUNCTION 


We consider the m = 1 case to represent the most ‘accurate’ model of volume exclusion effects 
for an on-lattice model of diffusion as it implies, simply, that one particle cannot overlap with 
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another. When coarse-graining this model, to consider random walk models with compartment 
capacities m > 1, our aim is that the mean and variance of particle numbers in each compartment 
are conserved. In what follows, we will sum the mean and variance of compartment occupancy of 
the accurate, m = 1, model over groups of m contiguous compartments to explore how accurate 
we can expect the coarse-grained model to be. 

To this end, we will define 


jm 


with the mean of and its variance. We wish to establish the relationship 

between: (i) and and (ii) v^^\t) and vj^\t). 

Steady state values. Eqs. Q-([^, together with the additional constraint that the sum of 
all variance and covariance terms must be zero (since N is constant), gives the steady states 
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where 1 < i, j < K and circumflexes are used to denote steady state values. It is then simple to 
check that 
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(15) 


i.e. the steady state means and variances are conserved through the process of coarse-graining. 


Time evolution. To obtain an expression for the evolution of // 
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We compare this to the coarse-grained model, Eq. (IM, which we re-state here for convenience: 
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To relate Eqs. Q and (16), and understand when we expect the coarse-grained model to replicate 
the dynamics of the accurate, m = 1, model, we need to establish a relationship between the 
(j = 1,..., iV) and the (j = 1, ■ ■ •, K). 

A natural choice for the coarse-graining would be to assume ~ /m for i = (j — 


l)m + 1,... ,jm. However, Eq. (16) would then give rise to a diffusion equation with constant 
diffusion coefficient m times larger than expected in the limit h —)• 0. This means that our coarse¬ 
grained model does not require the stringent condition that particles in a compartment of size m 
are uniformly distributed throughout that compartment. 

Instead, consistency between the accurate {m = 1) and the coarse-grained (m > 1) models 
arises from assuming particles in the m > 1 compartments are distributed, on average, according 
to a linear interpolation between m > 1 neighbouring compartments rather than being uniformly 
distributed throughout the compartment; 
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as 


shown in Eig. As a result we have 
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and evolution of the mean particle numbers in the coarse-grained system matches that of the 
accurate model. 



FIG. 2. Schematic of the interpolation process for to = 4. 


The entries of the covariances matrix cannot be interpolated in the same way, since its 

entries are positive on the diagonal and negative everywhere else. However, we can use similar rea¬ 
soning to argue that the variances and will also match, as presented in the Supplemental 
Material. 
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IV. NUMERICAL INVESTIGATIONS 


We now present numerical results to corroborate our findings. We consider the domain x G [0,1] 
with h = 1/128 and D = 10^. The initial condition is = 1, for j = 1,..., 16, and = 0 

otherwise, and attempts by particles to jump left out of compartment 1 or right out of compartment 
128 are aborted. 

We compare results from 5000 realisations of the random walk model with m = 1 with 5000 
realisations of the same model with m = 8 in Fig. The mean values predicted using both 
the m = 8 and PDF models are in excellent agreement with those predicted from the accurate 
m = 1 model. In addition, we see good agreement between the variances of the m = 1 and m = 8 
models |22] . Finally, we note that an additional advantage of the coarse-grained model is that 
generating realisations of the discrete random walk model with m > 1 can be achieved in 1/m^ 
of the time required by the m = 1 case, since the jump rates will be m? times smaller and so far 
fewer jumps will need to be simulated. 

To compute the error in the coarse-graining process, we solved the ODEs for both the means 
and variances, Eqs. Q and Q, over a range of values of m [23]. Results for m > 1 were compared 
against results with m = 1 using the histogram distance error (HDE) metric [24| : 

1 ^ 

HDE = - ^ |efc-pfcl, (19) 

^ A:=l 

where is the normalised value of the aggregated compartment of the m = 1 model and pk is 
the normalised value of the compartment of the m > 1 model. Figure [^demonstrates evolution 
of the error between the models with different values of m between t = 0 and t = 1. The HDE 
remains low in all cases observed, even though the initial condition does not satisfy the requirement 
that the densities in the m = 1 case can be interpolated onto the coarse {m > 1) lattice. 

V. CONCLUSIONS 

In this work we have used an on-lattice random walk approach to reconcile models of diffusive 
transport that incorporate the effects of excluded volume across spatial scales. Our work demon¬ 
strates that coarse-grained models incorporating simplified descriptions of excluded volume can 
be used in many circumstances, and these simplified models engender significantly lower compu¬ 
tational costs than their accurate counterparts. These computational savings are expected to be 
especially valuable for models in two or three spatial dimensions. However, care must be taken in 
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distance (x) 



distance (x) 

FIG. 3. Comparing the mean (a) and variance (b) of particle numbers for the m = 1 and m = 8 cases at 
t = 1, with K = 128/m. Dark grey bars: results from simulation of the random walk model with m = 1, 
light grey bars: results from numerical solution of Eq. (|^/Eq. Q with m = 8. Black dashed line in (a): 
solution of the limiting diffusion equation. 

pushing the coarse-graining process too far. For example, there is a delicate trade off between the 
initial conditions of the model and the size of m. Future work will be directed towards hybrid ap¬ 
proaches, in which a detailed description of the spatial dynamics can be retained where necessary, 
and the computational savings associated with the coarse-grained model taken advantage of where 
possible. 
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FIG. 4. HDEs for the means (a) and variances (b) in particle numbers as m is varied [22]. 

Appendix A: Derivation of mean and variance equations 

Recall that we have dehned a spatial lattice composed of K compartments. The distribution of 
particles over this domain is given by the vector n(t) = [ni(t), n2{t), ..., nxit)], where ni{t) denotes 
the number of particles in the compartment at time t. Each lattice compartment has capacity 
m, but in the interests of clarity we will omit the superscript (m) from the variables. Particles in 
compartment i may attempt to jump out to compartments i — 1 or i + 1: jumps in each direction 
are equally likely, and occur with rate per unit time of d. 

We define two operators, —>■ for i = 1,K — 1, and J~ : —>■ M^, for 

i = 2, ..., K, as 


4 

: [m,.. 

..rii, ....uk] 

[ni,.. 

— 2 ) 1 ? ^2 + l 1? ^2 + 2 ; •••; 
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: [ni,.. 
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Tli—2^ '^i — 1 I5 ^2 H“ 1? ^2 + 1 ; ^ 2 + 2 ; '^k]' 

(A2) 


Both operators move a particle into compartment i, taken from the compartment to the right or 
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left, respectively. We assume that attempted jumps into some compartment j fail with probability 
f{nj) due to exclusion effects. We can then write the probability master equation as 


d"’ ^ X] “ f{rii+i - 1)] Pr(J+ n,f) - n* [1 - /(rii+i)] Pr(n, t)} 

i=l 


K 


+ ^d{(ni + 1) [1 -/(ni_i - l)]Pr(J. n,t) - n* [1 -/(ni_i)] Pr(n,t)} . (A3) 
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Multiplying Eq. (A3| by nj and summing over all possible values that the vector n(t) can take we 
have 

M /K-l 

^ d{(ni + 1) [1 - /(ni+i - 1)] Pr(J+ n,t) - n* [1 - /(nj+i)] Pr(n, t)} 


^ = V n- 

dt ^ ^ , 

ni,n2,.--?^iC=0 \ ^=1 


+ '^d{{ni + 1) [1 - f{ni-i - 1)] Pr(J. n,t) - n* [1 - f{ni-i)] Pr(n,t)}^ .(A5) 
1=2 / 

We begin by considering only the first term of this expression, with Pr{J^ n, t) expanded explicitly 

to give 
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When i ^ j, j — 1 each term of this expression reduces to 


diujUi) - d{njnif{ni+i)), (A7) 


where {rijUi) indicates the mean of the product. When i = j, 
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Similarly for i = j — 1, 

Y dnj{nj-i + 1 ) [1 - f{nj - 1 )] Pr{ni,n2, ...,nj-i + l,nj - 1, ...nK,t), 

ni,n2,...,nj(=0 

= dl(nj-inj) - (nj-injf(nj)} + Mj_i - (nj-if(nj))]. (A 9 ) 
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We then consider the second term, 

N K-l 

- E E dui [1 - /(ni+i)] Pr{ni,n 2 , ■■■, rn, tikP), 

ni,n2,...,nK=0 i=l 

which evaluates to 

K-l 

- ^ d[{ninj) - {nif{ni+i)nj)] . 
i=\ 

When combined with the expressions derived from the first term these give us 

- Mj + (nj/(nj+i)) + Mj_i - {nj-if{nj)). 

Applying the same approach to the third and fourth terms, we obtain 
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Adding this expression to (A12), we arrive at 


dMj 

dt 


= d[-Mj + {njf{nj+i)) + Mj_i - {nj-if{nj))] 

+d [-Mj + {njf{nj-i)) + Mj+i - {nj+if{nj))] . 


(A14) 


Similar reasoning can be applied to derive equations for Mi and Mk and this approach can also 
be used to derive the variance and covariance equations, where we use 
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It can be shown that 

d ^ 
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By using {niUj) = Vij + MiMj, we can then write 
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with similar results obtainable for the covariance terms, and for i = 1 and i = K. 
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Appendix B: Matching variance equations 


We begin by noting that Eq. (8) of the main text is essentially a semi-discrete two-dimensional 
diffusion equation and has no dependence on m, except in determining the jump rate, so it is clear 
that, away from the diagonal of the covariance matrix, the covariances of the m > 1 model will 
evolve similarly to those of the summed m = 1 model. 

To consider the behaviour around the variance values, and adjacent covariances, we begin by 
noting that with m = 1, Eq. (7) of the main text becomes 

^ p [-4Vi + Mi_i (1 - Mi) + Mi (1 - Mi_i) 

+Mi(l-M+i) + Mi+i(l-Mi)], (Bl) 
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By the standard result for the variance of the sum of random variables, we find 


du 


(m) 


jm 


dt 


E 


dV 


( 1 )' 


dt 


jm-l jm /,..(!)' 

+ 2 E E 


p={j—l)m+l q=p-\-l 


dt 


(B3) 


We evaluate each of the terms on the right-hand side to obtain 
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We may draw an analogy to the mean equations, when after summation only those squares on the 
borders of the notional coarse grid compartment make a contribution to the dynamics; we wish to 
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FIG. 5. The red compartments represent the variances of the fully-excluding system (i.e. the central 
diagonal of the covariance matrix), while the blue compartments represent the covariances. It can be seen 
that particles are leaving the central compartment at every edge (accounting for the —term in the 
coarse grid model), and are entering from the surrounding compartments at each edge except for the two in 
the top-left corner and the two in the bottom-right corner (those adjacent to the fine grid variance terms). 
This explains why the input terms are — l)lm and (m — l)/rn, respectively. 


show that evolves in the same way as does in Eq. (7), 
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The contributions of the variance terms are illustrated in Figure It can be seen that they 
represent a simple diffusion pattern in two dimensions similar to the expressions seen in the mean 
equations. The only deviation from this is seen in the top-left and bottom-right corners, where 
there is outflow to the external cells but no inflow from them. Since the affected region on each 
edge is one compartment-length out of m, it seems intuitive that this would correspond to the 
m > 1 model, where the input terms from the surrounding compartments all scale with [m— l)/m. 
A similar argument can be made for the j terms, where outflow scales with (4 — 2jm). We 
also note that in our expression for Vj^\ the contributions from mean terms simply represent the 
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propensity for a move occuring, just as they do in the partially-excluding case. On this basis, and 
considering the strong agreements observed in Section IV of the main text, it seems reasonable 
to expect the variance of the total number of particles across m fully-excluding compartments to 
match that of a partially-excluding compartment with capacity m. 
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